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Abstract 

In the OPERA experiment, superluminal propagation of neutrinos can occur if one 
of the neutrino masses is extremely smaU. However the effect only has appreciable 
amplitude at energies of order this mass and thus has negligible overlap with the multi- 
GeV scale of the experiment. 



1 Introduction 



Recently the OPERA collaboration reported reported a measurement of the average time 
taken for neutrinos (f^ up to % level contamination) created at CERN (CN) to arrive at the 
Gran Sasso Laboratory (GS) compared to the time taken travelling at the speed of light in 
vacuo (c). They found an early arrival time of approximately 6t = 60 ns, which corresponds, 
at a significance of 6.0(7, to faster-than-light travel with speed {v — c)/c = 2.48 x 10^^ [1]. 

Since their announcement, a large number of papers have written that variously seek 
to explain it with or without new physics, question the experimental setup, or point out 
difficulties with new physics explanations (a selection is [21 [31 HI El E] ) . 

Here, we work entirely within the framework of standard relativistic quantum field 
theory and within the Standard Model (suitable modified to allow for neutrino masses and 
mixing matrix). We note that effective superluminal propagation of (both for individual 
events and, as we will later see, on average), is possible if the mass m of one of the mass- 
eigenstates is extremely small, so small in fact that space-like propagation can take place 
with appreciable probability even over the L = 730 km distance between CN and GS. 

To see that this is possible in principle, consider a neutrino of such a mass created at 
CN at time = and position x = 0, and arriving at GS with space-time coordinates 
{y^iV) = — 5t,L)^ To describe this we should use the Feynman propagator for spinors. 
However the neutrino beam has average energy (E) = 17 GeV [IJ. The neutrinos are created 
left handed (by weak decay of mesons in the decay tunnel) and being ultra-relativistic, will 
stay that way to very good approximation throughout their flight. Therefore we need only 
the left handed component, which effectively reduces the propagator to that of a scalar 
particle: 



(Here p'^ is the four-momentum.) When the interval s := {x^ — y^)^ — {x — y)"^ is negative 
as is measured by OPERA, i.e. is space-like, then we have that 



where Ki is a modified Bessel function. For large value of its argument it decays exponen- 

^We can neglect the two transverse spatial coordinates. From hereon we work in units with h — c = 1, 
translating back as necessary. 
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tially as: 



K^imV^s) ~ V-^ (3) 

^ (-s)i 



On the other hand for smah values, Ki diverges: 



Ki (4) 

vrmv— s 



We see that if we choose 

1 1 



m ~ 



3.8 X lO^^eV, (5) 



-s V2L6i 

then a substantial fraction of the neutrinos can propagate superluminally. 

Of course it is not the case here that the neutrinos are really tachyonic. The effect 
occurs because the neutrino remains off shellj^ The Feynman propagator arises from the 
Lorentz covariant time-ordered expectation of the and i>fj_ fields 

S{y-x) = {0\Tu{x)u{y)\0) 

= eix' - y'){0HxMym - eiy' - x'){o\u{yHx)\o). (6) 

The propagator arises as an intermediate step in a chain of reactions, for example being 
created at CN through the decay vr^ — t- fJ-^i^fj,, and absorbed at GS through f^n — t- /U^p. 

We can interpret this as the being created with energy-momentum p'^ = {E){l,v). 
The faster-than-light propagation occurs because the field operators are not localised but 
have a spread of order the Compton wavelength 1/m, which in our case is large enough to 
stretch from CN to GS. 

Since we know that — x^ > 0, i.e. that the neutrino arrives in GS after being created 
at CN, the first term currently plays no role. The potential problems with causality and 
Lorentz covariance arise when we view the propagation in a frame moving at speed vp 
sufficiently fast that the event at GS happens before the one at CN: y'^ = "ypiv^ — vpy) = 



7iry'^(l — vvp) < (where the dilation factor 'yp = 1/^1 — Vp). Since the Feynman 
propagator is Lorentz covariant, it still takes the same form, but in this frame it is the first 
term in ([S]) that now operates. It describes a super luminal antineutrino z^^ travelling from 
GS to CN. Energy and momentum conservation are preserved in the chain of reactions, and 



^It follows that the constraints noted in ref. |3J do not apply: kinematically electron-positron 
bremsstrahlung would here require momentum transfer with constituents of rock, and is just one of many 
diagrams heavily suppressed by the Fermi constant and the requirement that the neutrino remains off shell. 
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moreover we see that the neutrino energy has also reversed sign: E' = 'ypE(){l — vvp) < 0. 
Therefore we have a consistent interpretation in this new frame: a is created at GS with 
energy —E' and momentum pointing from GS to CN, through the process n — )• pfi^i>^, and 
absorbed at CN through the process P^7r+ — t- fi^. The kinematics of these processes are 
allowed because the anti-neutrino's 4-momentum is space-like. 

This is of course just the standard argument due to Feynman with small adaptations to 
this novel context, and demonstrates that the Feynman propagator allows an interpretation 
consistent with Lorentz covariance. However, there still remains an issue of causality. In- 
deed, in the new frame it seems the neutron at GS has somehow to know to decay beforehand 
so that the antineutrino can later be absorbed by the pion produced at CN. 

Once we have understood the process in more detail quantum mechanically, we can see 
how this issue with causality is resolved. We will return to this in the conclusions. 

While this paper was being prepared, ref. appeared, where similar ideas are proposed 
as an explanation for the OPERA measurement. In fact, if we adapt this observation to 
the setup in the OPERA experiment, we can see that the effect vanishes or at best is far 
too small to explain the measurement. 

The reason our setup is not yet the pertinent one is that it is not true that the neutrinos 
are created at an exact time. By Heisenberg's uncertainty principle, such a particle would 
have an infinite spread in energy. In sec. 2 we argue that the spread in energy for each 
neutrino's wave packet is 1/r ~ 0.2 eV about its central value Eq. In the ensemble of 
neutrinos that make up the neutrino bunches, these central values are spread over a wide 
range [H [7] but it is clear nevertheless that Eq is set at a much larger scale (GeV to multi- 
GeV). 

Given that each neutrino's energy is thus sharply peaked we can already expect that the 
effect must be heavily suppressed. From ([3]), the evanescent part of the propagator appears 
dominated by energies of order ([5]) - indeed we will see that it is a manifestation of quantum 
mechanical tunnelling, requiring energy to be less than m; this mismatch with the energies 
in the wavepacket ^ Eq ziz 1/t, ensures that any remaining effect is consigned to any small 
tail in the wavepacket that reaches down to these small energies. For example, if we assume 
that the wavepacket is a Gaussian of width 1/r then this supplies a suppression factor 

~ exp -T^Ei ~ exp -T^{Ef ~ (7) 

Alternatively, a wave packet with a lower cutoff > m on the neutrino energy would eliminate 
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the effect entirely. Furthermore, the space-hke neutrino then has to be detected at Gran 
Sasso, which requires overlap of ^ with the typical wavefunction for a detected neutrinojfl 
Such a wavefunction has again an order GeV central energy and spread ~ l/r. (At the 
most optimistic we can note that a lower bound is set by the lowest threshold energy for 
detection of z/^ via charge current interactions in the Gran Sasso Laboratory. This is for 
the process u^n — )• id~p, thus Eq > + ^m'j^/'mn = 112 MeV.) This results in multiplying 
by another suppression factor ~ exp —t'^Eq from the tail of its distribution. 

One might try to resuscitate the proposal by bringing l/r and Eq closer to each other. 
The relevant Eq could be as low as the threshold energy of 112 MeV and perhaps we have 
missed something in our estimate of the energy spread l/r in the neutrino wave packet 
and the true spread is far larger. If we for the moment take this optimistic view, we note 
that we still cannot fit the data. The dominance of the very low energy scale is not seen at 
OPERA: they repeated the analysis concentrating on only those charged current events 
occurring in the OPERA target (where reliable energies could be measured). They split 
this sample into bins of nearly equal statistics, taking events with energy higher or lower 
than 20 GeV. With a significance of greater than 3fT, they still see a super luminal velocity 
in the higher energy sample of the same magnitude (within errors) |1] . 

Another way one might try nevertheless to use this effect to explain the OPERA mea- 
surement is to boost ([1]). In other words, we note that the neutrinos are neither created 
at an exact time nor at an exact location; in reality we need to integrate over position 
space terms that supply the neutrino with the appropriate ultra relativistic momentum. 
One could then hope that the behaviour ([3]) would be the correct one for directions trans- 
verse to the neutrino's momentum, corresponding to small deficits in the energy required by 
the ultra relativistic on-shell energy-momentum relation. The amplitude would effectively 
take the form of the kilometres-wide wave function assumed in ref. [5], where the OPERA 
result is then explained without any adjustable parameters by the off-centre detection of 
these wave packets. However we will see that this set-up does not result in such transverse 
evanescence. 

Although neutrino species oscillate into one another, we can ignore this effect in the 
present paper. The known mass-squared differences [8] show that only one mass eigenstate 
can be as light as ([5]). Therefore the effect we are looking for is entirely due to evanescence 
^more accurately the amplitude constructed from the detected interaction products 
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in this mass eigenstate. The only effect of mixing is the multiplicative inclusion of mixing 
angles at the beginning where the CN converts to this state, and at the end where it 
reverts to f^, as detected in GS. 

The structure of the paper is as follows. In the next section (sec. 2) we estimate 
the initial spread in energy of the neutrino wavefunction at CN and the effective further 
decoherence due to its detection at GS. In sec. 3, we construct the initial wavefunction 
and the amplitude for neutrinos as seen at the Gran Sasso Laboratory, and draw out the 
piece relevant for this effect, confirming and extending the arguments given above. Finally, 
in sec. 4, we present our conclusions and also argue why these results ensure there is no 
violation of causality. 

2 Estimating coherence 

We have already noted that we cannot assume that the neutrino is created at a definite 
time. Likewise however, it is not true that each neutrino has been created with a definite 
energy. Such a particle would be completely delocalised in time, so the question of when it 
arrives in Gran Sasso would become meaningless. We therefore have no choice but to take 
into account the shape of its wavepacket, at least in its gross details, since the results will 
depend on this. 

It should be clear that here we are not discussing the energy spectrum of the ensemble 
of neutrinos in the beam, which is very broad, but rather the inherent quantum mechanical 
uncertainty in the momentum of a given neutrino when it is created at CN and any further 
inherent uncertainty imposed on it by being measured at GS. 

We follow closely the analysis in ref. [6j. We start at CERN with a proton bunch 
extracted from the Super Proton Synchotron. The energy uncertainty inherent in an indi- 
vidual proton wavefunction can certainly can be no smaller than that set by At ~ 5 ns, the 
smallest time features in the proton bunch [1]. This corresponds to Ai? = 1/ At w 1.3 x 10^'' 
eV. 

However, even if the proton beam is coherent at this level, it suffers decoherence on its 
way to becoming the beam. Firstly, the protons impact the graphite target, producing 
the mesons (mostly pions) that will decay to muon neutrinos. Initially these mesons are in 
a quantum state together with the other products of the collision (including various nuclei) , 
however they then suffer decoherence from thermalisation in a hot target, both directly and 
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also through their quantum mechanical coupling to the decay products. Assuming a target 
temperature of, say 300°C, this limits the energy-momentum resolution to ksT ~ 0.05 eV. 

Finally the mesons decay in the decay tunnel and here further decoherence takes place, 
again through coupling to the decay products (in this case the muon) . Consider for example 
the decay of a Tr"*" to ^i^Vp,. The resulting quantum state takes the form: 



where ^7r(q) is the wave function of the erstwhile pion, q and p are 3-momenta, and M 
stands for the matrix element for the decay. The muon is absorbed by a combination of rock, 
a Hadron stop and two muon detectors [1, 9j. This allows the experimenters to measure 
the transverse coordinates of the proton beam spot when it hits the target to a precision 
of ~ 50 - 90 ^m [S], however it is reasonable to assume that the rock itself localises the 
muons at the /xm level (similar to emulsion — see below), even if this is not recorded. This 
corresponds to a momentum decoherence of order 0.2 eV/c. Through ([5]) this decoherence 
is transferred to the neutrino. 

We conclude that the chief limiting factor on the coherence of the initial neutrino wave 
packet is through the decay of the mesons and results in a wave packet with an energy 
spread of l/r ~ 0.2 eV. 

For muon neutrinos that interact in the Gran Sasso detector, the impact spot is dis- 
cernible in the emulsion at the ^um level [10]; we can expect similar localisation in the rock 
in front of the detector for the external events. Therefore the act of measurement results 
in an effective energy spread in the wave packet of similar size to that in the initial packet. 

3 Tunnelling to Gran Sasso 

We can thus regard our neutrino as being created at position a; = at an uncertain time 
centred around t = 0, with energy localised to E = Eq with an accuracy l/r. Let us model 
the shape as a Gaussian. Then we have for the initial wavefunction: 




p phase space 




(8) 




(9) 



The amplitude to find the neutrino at Gran Sasso at time t2 = y is then 



oo 





(10) 



— oo 
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where S is the Feynman propagator ([T]). Since we are deahng with a situation where we 
know the neutrino arrives at Gran Sasso, we will choose the normalization factor so that 
the probability distribution reflects this: 

/oo 
dt2 \Af{t2) = 1. (11) 
-oo 

We can now make a further simplification. We can drop the transverse momentum inte- 
grations in ([1]) since these will only result in ~ 1/L^ losses that are scaled away when we 
normalise. Therefore, substituting ([1]) and using ([9]), we have 

exp {—iE[t2 — t]+ iqL} 



/oo roo 
dt "^{0,1) / dEdq 
-oo J —oo 



J^2 _ q2 _ _|_ 



exp{-iEt2-T^AEy2 + iqL} 
" li^^'^ E^-q^-m^ + ^e ' ^''^ 

where AE = E — Eq, and we have performed the t integration. Now we do the momentum 
integration. Since we know that L > 0, the ie prescription tells us to close the contour over 
the top. We obtain: 

-oo P 



where is the normalisation constant, where p = \J E'^ — rrfi has either a vanishing positive 
imaginary part (the ie) for \E\ > m, oi p = i\/m^ — E'^ for \E\ < m. 

We see that the exponential decay in ([3|) indeed arises from energies of order m, as we 
claimed. We also confirm that the exponential decay component is suppressed into the far 
tail of the probability distribution by ~ exp —t'^Eq. The amount of evanescent component 
thus depends crucially on the unknown shape of this tail. If we had chosen a wave packet 
with a sharp cutoff at some realistic minimum energy, we would eliminate the evanescent 
contribution completely. 

For completeness we note that (jl3p can be evaluated by the method of steepest descents. 
The dominant term comes from E ^ Eq giving: 

^(i2)|dommant ~ iV-^^exp- ft2 - — ) + i{Eot2 - PqL)] . (14) 

Pqt \2t^ \ vqJ 



(The approximation is valid providing \t2 — L/vo\ ^ t^Po ) Here pq = ^ Eq — m? is the 
central momentum and vq = pq/Eq is the classical velocity. We recognise this as nothing 
but the expected result of propagating the wave packet to Gran Sasso without dispersion. 
The evanescent part can also be evaluated. Writing E = m(l — z) for small positive z at the 
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top boundary, the integral can again be evaluated by steepest descents. Adding to this the 
term from the bottom boundary E = —m{l — z), and similar size pieces from E = ibm(l + z) 
(which can again be evaluated by steepest descents), we get for the evanescent part 



M ,T / 27r _^2j72/9 \ niLp' vr 1 , , 

A{t2)\evanescent » -2lN J —6 ^ ^o/^gxpz <^ - mt2 + - > . (15) 

This approximation is valid providing mt2 ^ 1- We see again the damping by the tail of 
the wave function. We can now carry this through to a computation of the superluminal 
component of velocity. The evanescent term will provide one, but we do not present the 
computation since we have seen that it necessarily depends on a vanishingly small unknown 
quantity. 

Finally, we address the question raised in the introduction: whether, after taking into 
account appropriate spatial dependence of ^ and spatial integration in (jlOp . so as to incorpo- 
rate the fact that the neutrino has some momentum po (with associated small uncertainty) 



which is slightly larger in magnitude than the energetically allowed po = \J Eq — m?^ the 
resulting amplitude could take the form similar to ()14p in the direction of po but with trans- 
verse evanescent tails corresponding to imaginary transverse momenta? Recall that such an 
amplitude would be very similar to the wave packet envisaged in ref. [5] where results con- 
sistent with the OPERA measurement were derived as a result of off-centre measurement 
of these wave packets. 

Clearly in order to investigate this we now need to keep the transverse momentum 
integrations in ([TJ . The modified wave function and spatial integration incorporated in (jlOp 
results in a momentum integral in ()12p of the form: 

,7(L) = (16) 

J i/^ — — -|- le 

where strongly peaked about q = po incorporates the momentum dependence induced 
by the initial wave function. 

Now we appeal to the Grimus-Stockinger theorem jllj which states thal|f 

J(L) = (-PL/L) exp{ipL} + 0(^-^/2), (17) 



if p = y/W — m? is real, and is O^i/L?) otherwise. We see that providing E > m we 
get only the analogous dependence to ([H]). Once we integrate over energy in ([T2]), the 



^providing $ satisfies the reasonable conditions that it is three-times continuously differentiable and first 
and second derivatives decrease at least as for |q| — !> oo 
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mismatch between pQ and pQ will result in suppression arising from the partial overlap of 
exp — r^Ai?^/2 and <I> {—pL/L), however we do not generate evanescent tails as a result of 
this mismatch. On the other hand, for E < ni, J decays faster in L, which we can associate 
with the evanescence. We see that even with this 'boosted' wave packet, it is still the case 
that the evanescent behaviour responsible for superluminal propagation, is restricted to the 
regime E < m. 

4 Conclusions 

We have seen that although superluminal propagation is possible over the 730 km from 
CERN to Gran Sasso if the mass of the lowest neutrino mass eigenstate is so small that it 
remains off shell, cf. dS]), the effect does not survive projection on the relevant energy scales, 
being killed dead or by typically a huge suppression ([7D. This projection is required because 
one must integrate over the initial neutrino wave function which carries a very rapidly 
oscillating exponential, set by the multi-GeV energy of the created neutrino. One could 
reduce the suppression [i.e. increase the overlap with energy scale m) by concentrating on 
muon neutrinos with energies only just higher than the threshold energy of 112 MeV. This is 
still not enough: one would have to argue that the quantum-mechanical energy uncertainty 
in the neutrino wave packet was not of order 0.2 eV as we have argued but - implausibly - 
tens or hundreds of MeV. Even if one does this, the effect is still concentrated at energies 
of order m, resulting in the wrong energy dependence. 

One could reduce the time for which the neutrino has to remain off shell by considering 
cases where it propagates on shell to a point near the detector and, due to interaction with 
the rock, is then kicked off shell for the remaining part of its journey. Note that an off-shell 
neutrino can be passed from this point to the detector in principle almost instantaneously. 
In this way we achieve the smallest negative interval s from the instantaneous jump (0, 5t) 
from {L — 6t,L — 6t) to the detector. In ([5]) we then get the larger mass m ~ l/5t ~ 10~^ eV 
[3]. This is still far too small to bridge the gap to the energy scales set by the experiment. 

Although the OPERA measurements cannot be explained by assuming very low mass 
off-shell neutrinos, it is still the case that such neutrinos can propagate superluminally. How 
can this be reconciled with causality, in particular why does this not lead to faster-than- 
light communication? Note that we have shown that these neutrinos must carry energy E 
less than m. Leaving aside practical issues involved in detecting such neutrinos, we note 
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that even in principle, a detector for such neutrinos would have to be restricted to using 
wavelengths larger than or of order their Compton wavelength A = otherwise the 

observation itself would disturb the system too much - for example by pair creating the 
very neutrinos it was trying to measure. Thus we see that the restriction to E < m for 
super luminal neutrinos, which comes out from the detailed analysis, is in fact necessary to 
ensure that they cannot be detected with sufficient spatial resolution to allow faster than 
light signals. 
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